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Abstract 

Vortex structure of pure d^2_y2-wave superconductors is microscopically an- 
alyzed in the framework of the quasi-classical Eilenberger equations. Self- 
consistent solution for the d-wave pair potential is obtained for the first time 
in the case of an isolated vortex. The vortex core structure, i.e., the pair 
potential, the supercurrent and the magnetic field, is found to be fourfold 
symmetric even in the case that the mixing of s-wave component is absent. 
The detailed temperature dependences of these quantities are calculated. The 
fourfold symmetry becomes clear when temperature is decreased. The local 
density of states is calculated for the selfconsistently obtained pair potential. 
From the results, we discuss the flow trajectory of the quasiparticles around 
a vortex, which is characteristic in the d^2_y2-wsive superconductors. The 
experimental relevance of our results to high temperature superconductors is 
also given. 

PACS numbers: 74.60.Ec, 74.72.-h 
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I. INTRODUCTION 



A number of investigations were carried out theoretically and experimentally to iden- 
tify the symmetry of pairing state in high-T^ superconductors. Although precise pairing 
symmetry has not been determined yet, it is recognized that (i^2_y2-wave symmetry is most 
probable.0 Recently, the vortex structure of the d-wave superconductors attracts much at- 
tention because it may have different structure from that of conventional s-wave supercon- 
ductors. The internal degrees of freedom in d^-i-yi-wave pairing, i.e., k"^ — ky in reciprocal 
space, has fourfold symmetry. It is expected that reflecting the symmetry, the core structure 
of isolated vortex may break the circular symmetry and show fourfold symmetry in d^^-y^- 
wave superconductors. In conventional s-wave superconductors, the winding number of the 
pair potential is 1 around the vortex. In (i^2_y2-wave superconductors, other components of 
winding number 4n + 1 (n: integer) may be induced and make the vortex structure fourfold 
symmetry. 

Fourfold symmetric vortex in dx2_y2-wave superconductors was so far mainly considered 
by Ren et al.ii and Berlinsky et al.i on the two-component Ginzburg-Landau (GL) theory 
for s- and rf-wave superconductivity. According to the consideration based on the two- 
component GL theory, it is possible that the s-wave component is coupled with the d-wave 
component through the gradient terms. Therefore, the s-wave component may be induced 
when the d-wave order parameter spatially varies, such as near the vortex or interface under 
certain restricted conditions.!'! The induced s-wave component around the vortex is fourfold 
symmetric. The resulting vortex structure in d-wave superconductors, therefore, exhibits 
fourfold symmetry. The structure of the induced s-wave component was also analyzed by 
solving the tight-binding Bogoliubov-de Gennes equation self consistently on a 16 x 16 latticell 



and by the quasi-classical Eilenberger theoryJH 

However, this fourfold anisotropy clearly appears only when the amplitude of the induced 
s-wave component is comparable to the (i^2_j^2-wave component. It is still not clear both 
theoretically and experimentally whether or not this is indeed the case of high temperature 
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superconductors. In the case where the s-wave pairing interaction is neghgibly small com- 
pared with the (i^2_j^2-wave one, the induced s-wave order parameter is negligible. In the 
limit As 0, the GL equation reduces to an equation for only (i^2_j^2-wave component, 

— + dy)A,j.2^y2 — A^2„j^2 + |A^,2_j^2pA^2„j^2 = (1.1) 

in the usual dimensionless form. Since Eq. (|l.l| ) is same as that of the conventional s-wave 
superconductors, the vortex structure remains circular symmetric within the GL framework. 

Strictly speaking, the GL equation (|1.1|) is valid only near the transition temperature. 
Far from the transition temperature, we have to include several correction terms. Among 
them, there are some terms leading to fourfold symmetry, such as the so-called non-local 
correction term,i'i 

- 7ln (^) [lidl + dlf - a^aj}A.2_,2, (1.2) 

where 7 = 62C(5)/49C(3)^. Including these terms neglected in the conventional GL theory, 
we can investigate the fourfold symmetric vortex in pure rf-wave superconductors. On low- 
ering temperature, these correction terms make important role and fourfold symmetry is 
expected to be evident. 

The purpose of this paper is to investigate the vortex structure in pure c?2.2_j^2-wave 
superconductors using the quasi-classical Eilenberger equations,! and show that fourfold 
symmetric vortex is realized even in the case the induced s-wave component is absent, i.e., 
in the case that the s-wave pairing interaction can be neglected. Merits of the quasi-classical 
calculation are as follows. It can be applied at arbitrary temperatures, while the GL theory 
is valid only near the transition temperature. We can consider the terms neglected in the 
GL theory, the higher order terms of order parameter and the higher order derivative terms. 

The quasi-classical calculations on the vortex structure were selfconsistently carried out 
for conventional s-wave superconductors by Pesch and Kramer,0 and Klein0 quite thor- 
oughly. For pure (i2.2_j^2-wave superconductors, the quasi-classical calculations on the vortex 
were so far performed by assuming circular symmetry for the amplitude of the c?2.2_j^2-wave 

3 



pair potential. Schopohl and Maki showed that local density of states exhibits characteristic 
fourfold symmetric distribution.lll The current authors demonstrated that the supercurrent 
and magnetic field distribution around an isolated vortex are fourfold symmetrical These 
results suggest that the pair potential itself is fourfold symmetric although the circular sym- 
metry is assumed in the initial pair potential. In this paper, the selfconsistent calculation 
for the (i2;2_j^2-wave pair potential is performed for the first time. The fourfold symmetry 
of vortex structure, i.e., the pair potential, the supercurrent and the magnetic field distri- 
bution around a vortex is demonstrated to exist and the temperature dependence of these 
quantities is investigated. We also calculate the local density of states in detail based on the 
selfconsistently obtained pair potential. 

Here we consider the case of an isolated vortex under a magnetic field applied parallel to 
the c-axis (or z-axis) in the clean limit. The Fermi surface is assumed to be two-dimensional, 
which is appropriate to high- Tc superconductors, and isotropic for simplicity. Throughout 
the paper, energies and lengths are measured in units of the uniform gap Aq at T = and 
the coherence length = vp/ {vp'- Fermi velocity), respectively. 

The rest of this paper is organized as follows. In Sec. we describe the method of 
calculation based on the quasi-classical Eilenberger theory. Section |in| contains selfconsistent 
numerical results of the vortex structure computed at various temperatures. The summary 



and discussions are given in Sec. IV 



II. QUASI-CLASSICAL EILENBERGER THEORY 

To obtain the quasi-classical Green functions, we solve the quasi-classical Eilenberger 
equations for a pair potential, 

A(^,r) = A(e,r)e^^ (2.1) 

where r = |r| = y/x'^ + y'^ is the distance from the center of a vortex line situated at 
the origin, 6 is the angle of k- vector with the a-axis (or x-axis), and e*"^ = {x + iy)/r is 



4 



factored out for later convenience. After an appropriate gauge transformation about 0, the 
Eilenberger equations for the quasi-classical Green functions with the Matsubara frequency 
ujn = (2n + 1)'kT are given as 

{^n + \ {d\\ + ^^110) }/>„, r) = A(e, r)(^K, 9, r), (2.2) 
[un - \{d\\ - id\\ct)) ]f\un, 9, r) = A*(^, r)^K, 9, r), (2.3) 
d\ig{u;^, 9, r) = A*{9, r)/>„, 9, r) - A{9, r)/t(^„, r), (2.4) 

g{uJn, 9, r) = (l - ^, t)P{uj^^, 9, r)) Re(7(cu„, 9, r) > 0, (2.5) 

where d\\ = d/dr\\ and d\\(l) = —r±/r'^. Here, we have taken the coordinate system: u = 
cos ^^x + sin 9'y, v = — sin 9^ + cos ^^y, thus a point r = sic + is denoted as r = ryu + r_LV. 
The anomalous Green functions / and in Eqs. (p.2|) - (|2.5|) are related to the usual notations 
/ and as / = fe^'^ and = /"''e"*'^. Since we consider the isolated vortex in the extreme 
type II superconductors, the vector potential in Eqs.( |2.2| ) and ( |2.3| ) can be neglected. 

We solve the first-derivative equations ( |2.2| )-( p^ ) along the trajectory where r± is held 
constant. To obtain the quasi-classical Green functions. 
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(2.6) 



we use the so-called explosion method. In addition to a physical solution ^p^, Eqs. 
(p.2|)-( pl4D have two unphysical solutions ^+ and The solutions g± explode (increase 



exponentially) in the directions ±k and decrease in the opposite directions. Even when we 
use the physical solution as a initial value, the unphysical solutions always mix and become 
dominant during the process of the numerical integration of Eqs. p.2| )-( [2l^ ) along a long 
path. We obtain g± by integrating from ry =prA to ry, where rA(> 0) is large so that explosion 
takes place. It is known00 that the physical solution is obtained from the commutator of 
the two unphysical solutions, 



9ph = c[g+,g-], (2.7) 

where c is a constant determined from Eq. (|2.5|) . 

The quasi-classical Green functions are also obtained by the method of the Riccati trans- 
formation. The transformation removes the unphysical solutions in the numerical integra- 
tion. Using the parameterization devised by Schopohl and Maki,@ 

J 2a p 2b 1 — db 

1 + db^ 1 + db^ 1 + a6 ' 

in Eqs. (|2.2|)-(P^, we obtain the Riccati equations 

d\\d{un, e, r) = A{e, r) - {2ujn + + A*(^, r)a(u;„, 6, r)}a(cu„, 6, r), (2.9) 

9||6(^„, e, r) = -A*(^^, r) + {2^„ + + A{e, r)6(^„, 9, r)}6(^„, 9, r). (2.10) 
We integrate Eqs. ( |2.9| ) and ( p.lO| ) using the solution far from the vortex 



iol + \Ai9,rW-uj^ ^ujl + \Ai9,vW-uj^ 

boo = T77r~\ ' (^n > 0) (2-11) 



A*{9,r) ' A(^^,r) 

as an initial value, and obtain a, b and the quasi-classical Green functions. We confirm that 
both methods, the explosion method and the Riccati transformation method, give the same 
solution. 

The selfconsistent condition is given as 

f 2^ d9 . 

where Nq is the density of states at the Fermi surface. The pair potential and the pairing 
interaction are decomposed into s-, dx2-y2- and ci^^y-wave components, 

V{9,9') = Vs + V^2_y2cos{29)cos{29') + V^ysm{29)sm{29'), (2.13) 
A{9, r) = A,(r) + A^,2_j^2(r) cos(2^) + A^,j,(r) sin(2^). (2.14) 



A(^,r)=iVo27rT 5: / —V{9,9')Ruj^,9\v), (2.12) 



Substituting Eqs.( |2.l3| ) and ( f^.l4| ) into Eq.(^TT2|), we obtain the following selfconsistent 
equations for each component: 
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A,(r) = VsNo2nT ^ / — /K, 6, r), (2.15) 

JO iTT 



a;„>0 



/O 277 



A,2_y2{T) = V,2_y2No2nT / — />„, ^, r) cos(2^), (2.16) 



/•27r _ 

A,,(r) = V,yNo2nT J2 / 7r/(^"' ^' sin(2^). (2.17) 



U)n>0 



In Eqs. (|2.15| ) and ( p.l7| ), As(r) and Axy{T) are proportional to Vg and V^^, respectively. 



Here, we consider the case \ Vs\, \ Vxy \ -C V^2_y2 as mentioned before. In this case, since As(r) 
and Axy{r) are negligible compared with A^2_y2{r), the s- and c?^.y-wave components do not 
affect (ij,.2_j^2-wave superconductivity. For the case Vs and V^y are negligible but finite, the 
structure of As(r) and A^y{r) are calculated from Eqs. ( p.l5| ) and ( |2.17|) respectively, which 
was in detail reported in our previous paper .0 



From now on, we set As{r) = A,j;y{r) = 0. Equation ( |2.16| ) is the selfconsistent condition 



for (i^2„y2-wave superconductivity. In our calculation, we use the relation 

^ log^ + 27rT A' (2.18) 



V^2_y2No Tc o<uj^<^c 1^' 



and set the energy cutoff: ujc = 20Tc. We calculate the r.h.s. of Eq. (|2.16| ) and obtain 



the new value for Aj.2„j,2(r). Using the renewed pair potential, we solve the Eilenberger 
equations (p.2|) - (|2.5|) again. Starting from the initial form 

A^2_j^2(r) = A(T) tanhr, (2.19) 

we repeat this simple iteration procedure twenty times and obtain a sufficiently selfconsis- 
tent solution for A^2_y2{r). Here, A(T), the temperature-dependent uniform gap without 
magnetic field, is obtained from the BCS relation. 

In our numerical calculations, we discretize r in the region r < rc = 10 for Aa,2_j^2(r). 
In order to avoid the spurious symmetry breaking due to computational artifacts, the mesh 
points are located on the cylindrical coordinate so that the choice of mesh points is designed 
not to break cylindrical symmetry. When we solve Eqs. (|2.2|) - (|2.4| ), we need to know 



Aa;2_j^2(r) for arbitrary r. It is given by interpolation of the values on the mesh points. Far 
from the vortex r > r^, we assume the pair potential to be the form of Eq. ( |2.19| ). To confirm 
the validity of this assumption, we calculate also for the other rc or other asymptotic form, 
such as A(T){1 — (2r^)~^}. These changes do not affect the final results of the vortex core 
structure. 

The supercurrent around a vortex is given in terms of g{uJn, 6', r) by 

J(r) = 2et;i.iVo27rT ^ / ^^(cu„,^,r). (2.20) 
The associated magnetic field distribution is determined through 

— J(r) = V X H(r). (2.21) 

c 

The local density of states is given by 

N{r,E)= Reg{iu;^^E + zr],9,r), (2.22) 

JO 2,71 

where 77 is a positive infinitesimal constant. To obtain g{iujn E + iri,6,r), we solve 
Eqs.( |2.2| )-(|2.5|) for r] — iE instead of u;„ using the selfconsistently obtained pair potential. 



III. FOURFOLD SYMMETRIC VORTEX STRUCTURE 

A. Pair potential 

We calculate the vortex structure at T/Tf. = 0.1 ~ 0.7. First, we consider the dx2_y2-wave 
pair potential. The pair potential A^2_y2[r) at T/Tc = 0.1 is shown in Fig. |l]. The amplitude 
of the pair potential clearly exhibits fourfold symmetry as shown in Fig. |1] (a). The amplitude 
along the 0° direction (along the a;-axis and y-axis) is suppressed compared with that along 
the 45° direction (along the lines y = ±x). The phase of the pair potential also exhibits 
fourfold symmetry, as shown in Fig. |1| (b). In conventional s-wave superconductors, since 
A is real, argA = arg(Ae*'^) = 0, that is, the phase increases uniformly around a vortex. 
In contrast, the phase does not increase uniformly in (ij.2_j^2-wave superconductors. The 
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phase argA is larger than in the region < < 7r/4, and smaller than in the region 
7r/4 < < 7r/2. These behaviors of the pair potential can be explained as follows near the 
vortex center, where the pair potential can be expanded in the odd power of r. For the case 
of the fourfold symmetry, the pair potential consists of terms with e*^^""'"^)'^ (n: integer). 
Following the well-known consideration about vortex core,lii the small r-expansion begins 
from the order r^'^"'+^\ for the factor of e**^^""''^)'^. The pair potential is, therefore, given as 
follows near the vortex center, 

A,2_y. = A,2_y.e'^ = ao{(r + bor^)e''^ + Cor'e-''^ + ©(r^)}, (3.1) 



where ao, &o and Cq are constants. From Eq. ( |3.1| ), the amplitude and phase of the pair 
potential are, respectively, obtained as follows, 

\A^2_y2\ = ao{r + bor^ + cor^ COS 4:(f) + 0{r^)}, (3.2) 

argAa.2_y2 = — cor^ sin40 + ©(r"^). (3.3) 

The case Cq < is consistent with our numerical calculations shown in Fig. Far from the 
vortex, the fourfold symmetry gradually fades and reduces to circular symmetry. 

The temperature dependence of Aa.2_y2(r) is presented in Fig. |^. Figure ^ (a) is for the 
amplitude of the pair potential as a function of r along the 0° direction, |A^2„j^2(r, = 0)|. 
It shows that the size of the vortex core steeply decreases on lowering temperature. Also in 
the conventional s-wave case, this reduction was reported by Pesch and Kramer ,0 and Gygi 
and Schliiter.lli Concerning to the temperature dependence of fourfold symmetry. Fig. |^ 
(b) shows the difference of the amplitude along the 0° and the 45° direction, |A^2_y2(r, = 
0)1 — I A^2_y2(r, = 7r/4)|. And Fig. ^ (c) shows the phase of A^2_y2 along the 22.5° direction, 
arg A^.2 _j^2 (r, = it/ 8). From both figures, we see that fourfold symmetry of vortex becomes 
clear on lowering temperature. To examine these temperature dependences quantitatively, 
the curves in Figs. (a) and (b) are fitted by the following functions respectively in the 
core region r < tq/S, where tq is defined as the radius r giving minimum in Fig. |^ (b) for 
each curves, 

9 



/\3 /\5 /\7 

ly / 7^ \ (fx (if 



|A,2_,2(r,0 = O)|-|A,2_,2(r,0 = 7r/4)|=A(T)<j6i- + 63(7) +^5(7) +^7(7 



(3.5) 



where ^, as, as, ay, 61, 63, 65 and 67 are fitting parameters. We attain the best fit with 
the values summarized in Table |. The obtained value for hi is very small and regarded as 
0, which is consistent with the consideration leading to Eq. ( p.2|) . The fitting parameter 
^ thus obtained characterizes the coherence length which corresponds to the actual size of 
the vortex core. Its temperature dependence is plotted in Fig. |^ (a), showing that the core 
size ^ linearly approaches on lowering temperature. From Table |, we can see that the 
contribution of the higher order terms increases with decreasing temperature. Figure ^ (b) 
shows the temperature dependence of the fitting parameters aa and 63. When temperature 
is decreased, both |a3| and I63I similarly increase, which means that the contribution of 
the fourfold symmetric term increases. On the other hand, in the limit T Tc, the 
parameters approach the values expected from the GL equation, ~(8'Cgl)^^5 ^3 ~^ 0- 

B. Supercurrent and magnetic field around a vortex 

We show a stereographic view of the current distribution in Fig. |^ (a). It is clear that the 
amplitude of the current |J(r)| is fourfold symmetric. There are four small peaks around 
the core at 45° and its equivalent directions. The breaking of the cylindrical symmetry 
is stronger in the core region, while toward the outer region the cylindrical symmetry is 
gradually recovered. To consider the direction of the current flow, we decompose the current 
to the radial and the rotational components. The rotational component has almost the 
same distribution as |J(r)| shown in Fig. ^ (a). The radial component Jr(r) is shown in 
Fig. ^ (b), which indicates that the current has small fourfold symmetric radial component. 
The counter-clockwise current has negative (positive) Jr and curves inward (outward) for 
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O<0<7r/4(7r/4<0<7r/2). The current flow, therefore, deviates from circular trajectory, 
and forms fourfold symmetric trajectory, which comes near the center of vortex at = 7r/4. 

The associated magnetic field distribution H{r) is displayed in Fig. ^. The contours of 
constant magnetic field coincide with the supercurrent streamlines. Reflecting the current 
distribution, the magnetic field is also fourfold symmetric. In the core region, the field 
extends along the 0° direction, rotated by 45° from the current distribution. This is because 
the field is strongly screened by the induced current along the 45° direction. 

Next, we consider the temperature dependences of these quantities. The current and the 
magnetic field distribution are shown, respectively, in Figs. ^ and at various temperatures. 
When temperature is lowered, the peak of the current moves toward the core and its height 
increases, as shown in Fig. |^ (a), and magnetic field distribution shrinks toward the core, 
as shown in Fig. ^ (a). These features reflect the narrowing of the core size with decreasing 
temperature as mentioned before. As for the anisotropy, the fourfold symmetry becomes 
clear in the current, as shown in Figs. ^ (b) and (c), and the magnetic field, as shown in Fig. |^ 
(b), with decreasing temperature. It is seen from Fig. |^ (b) that | J(r, = 0)| — |J(r, = 7r/4)| 
changes sign, positive outside of the vortex core. Because of this change, the fourfold 
symmetric magnetic field becomes almost circular beyond the core region. 

C. Local density of states 

The local density of states A^(r, E) is calculated for the selfconsistently obtained pair 
potential. We confirm that the selfconsistent calculation gives essentially the same fourfold 
symmetric local density of states as that predicted by Schopohl and Makilll using a test 
potential A^2_y2(r) = A(T)tanhr. Figure ^ shows N{r,E) for E=0.2, 0.01 and where 
we use the selfconsistent pair potential at T/Tc = 0.1. From Fig. ^(a) it is recognized that 
A^(r, E) for E = 0.2 has the similar fourfold symmetric distribution to that of Schopohl and 
Maki in the core region. Because of narrowing of the core size at low temperatures in our 
selfconsistent calculation, the characteristic length scale for the variation of A^(r, E) is much 
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shorter than that of theirs. Figure ^(b) shows A^(r, E) for E=0.2 in a wider region. It is 
noted that the eight ridges of A^(r, E) extend to far from the vortex. This feature is not 
covered by Schopohl and Maki.0 To consider the case of lower E and the hmit E ^ 0, we 
show iV(r, E) for E = 0.01 in Fig. |(c), and for E = in Fig. |(d). 

The distribution of A^(r, E) can be understood as follows. The peak lines in Fig. |^ 
consist of four flow lines of quasiparticle with energy E. They screen the magnetic field in 
total. The flow lines are schematically presented in Fig. where four lines are labeled as 
1~4, respectively. The quasiparticle on the line 1 starts toward the 45° direction (the node 
direction of the dj.2_y2-w8ive superconductivity) at far from the vortex, and approaches the 
vortex. It turns around at the vortex core in a parabolic orbit, and departs from the vortex 
core toward the node direction of —45°. In this flow trajectory, the flow direction changes 
from 45° to —45° gradually so that the factor cos2^^ of the (i-wave pair potential does not 
change its sign. The four flow lines of quasiparticles around a vortex are the characteristic 
point of the (i^2_„2-wave superconductivity, which is contrasted with the circular flow§ m 
the conventional s-wave superconductivity. From Fig. |^ (b)-(d), we recognize that the flow 
trajectories in Fig. ^ approach the center of the vortex with decreasing E, and reduce to the 
lines y = ±x in the limit E 0. The A^(r, E) for E = in Fig. ^ (d), therefore, has four 
small ridges in the 45° and its equivalent directions, in addition to the usual large zero bias 
peak at the vortex center. On the other hand, when E becomes larger, these small ridges in 
the outer regions are invisible because the quasiparticles with larger E are in the scattering 
state {E > Aocos26'), thus they distribute uniformly. 

In Fig. |10| (a), (b) and (c), we show N{r,E) as a function of E and distance from the 
vortex center along the directions = 7r/8, = and = 7r/4 respectively. This type of 
figures is already calculated by Schopohl and Maki0 who evaluated it non-selfconsistently. 
We also note that Hess et al. and Renner et al.S performed the STM experiments to 
directly observe the local density of states A^(r, E) in a s-wave superconductor 2H-NbSe2 
whose characteristic features are analyzed theoretically by Gygi and Schliiter,0 and Shore 
et al.lii Three of the four flow trajectories in Fig. ^ can be seen as three peaks for E < 1 in 
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Fig. [ly (a). To specify their relation, we label the three peak lines and the corresponding 
positions as A, B and C in Figs. ^ and ^ respectively. For A^(r, E) along the direction 
= 0, the peak lines A and B overlap each other as seen from Fig. |10| (b). With increasing 
(f) from to 7r/4, peak lines A and C shift to lower energy and line B shifts to higher energy. 
It is seen from Fig. |TU| (c) that for N{r,E) along the direction = 7r/4, the peak lines B 
and C overlap each other and the line A reduces to the peak at E = 0. The peak line at 
E = 1 in Fig. |l^ corresponds to the gap edge, which becomes clear far from the vortex. 

IV. SUMMARY AND DISCUSSIONS 

The vortex structure in a pure dx2-y2-v/ave superconductor is studied in the framework of 
the quasi-classical Eilenberger equation. By a selfconsistent calculation, we obtain the pair 
potential around an isolated vortex, which shows the characteristic fourfold symmetry in the 
core region. The associated supercurrent and magnetic field also exhibit fourfold symmetry 
around a vortex. We confirm that the fourfold symmetry becomes clear with decreasing 
temperature. Using the selfconsistently obtained pair potential, the fourfold symmetric 
local density of states is calculated. 

While our numerically obtained vortex structure has similar fourfold symmetry to that 
predicted by Xu et all who use the two-component GL equations, the origin of the fourfold 
symmetry is quite different from theirs. In their theory, the fourfold symmetry is induced by 
the mixing of s-wave component. On the other hand, we consider the case where the mixing 
of the s-wave component is negligible, i.e., the pure (i^2_j^2-wave case. Even in this case, the 
vortex structure exhibits fourfold symmetry. It is due to the fourfold symmetric terms which 
are higher order of ln(T/Tc) and neglected in the conventional GL theory. When temperature 
decreases below Tc, the contribution of these terms increases and fourfold symmetry becomes 
clear. Since our quasi-classical calculation automatically takes these higher order terms into 
account, we can obtain the fourfold symmetry of the vortex structure even in the pure 
dx2_y2-vf&Ye superconductors. 
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The scenario based on the two-component GL theory mentioned before is not reahzed 
in the case the induced s-wave component is neghgible. In our formulation, it corresponds 
to the case where the s- (or extended s-) wave pairing interaction Vg is neghgible compared 
with Vx2-y2 in Eq. ( |2.13| ). The induced s-wave component As(r) is proportional to Vg in 
Eq. ( |2.15| ). In such a case |K| ^ Vx'^-y'^, ^s(r) is negligible and does not effect the dr^2_y2- 
wave superconductivity. In this case, our scenario based on the quasi-classical theory is 
realized. The consideration of the limit 14 on the GL theory gives the same result, if 
the Fade approximationii is not used in the derivation of the GL equations. The coefficient 
of A<j in the GL equations is given by = {VsNo)~^ — 2{Vx2_y2No)^^ without the Fade 
approximation. Since \As\ oc A^ ^ in the limit Vg 0, which is consistent 

with our consideration. On the other hand, in the GL theory with the Fade approximation, 
tts = 2(V^2_y2Ai"o)^^{l + 2(— 14)/V^2_j^2}, which gives an invalid result ag finite in the limit 
Vg ^ 0. The Fade approximation is inadequate in the case V^^ — and Vg is repulsive. When 
Vg is repulsive, the original GL theory without the Fade approximation gives unphysical 
solution for A^.ii To understand the repulsive case, further careful studies are needed. A 
selfconsistent quasi-classical calculation may be one of the possible approaches. 

The four flow trajectories of the quasiparticle for E{< 1), as shown in Fig. ^ schemat- 
ically, is the key point to understand the fourfold symmetry of the vortex structure in a 
dj;2_y2-w8ive superconductor. It is contrasted with the circular flow in conventional s-wave 
superconductors. The local density of states in Figs. || and |10| directly reflects the four 
flow trajectories. As the supercurrent is the sum of the quasiparticle flow contributions with 
various £"s, the resulting supercurrent and the associated magnetic field distribution around 
a vortex have fourfold symmetry. Fourfold symmetry of the pair potential also reflects the 
quasiparticle distribution. 

Reflecting the fourfold symmetry of the vortex core, the vortex lattice or flux line lattice 
may be distorted from hexagonal symmetry in ci-wave superconductors. The distortion 
may occur at higher field as the inter-vortex distance is short and the effect of the fourfold 
symmetric core structure increases. In the case of anisotropic core structure, we have to 
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consider the orientation of the vortex lattice, i.e., nearest neighbor direction of vortices 
relative to the underlying crystal lattice. Won and MakE considered that at near Hc2 
a square lattice tilted by 45° from the a-axis is preferable to a triangular lattice at T < 
0.8Tc in (i-wave superconductors. By using the two-component GL theory mentioned above, 
Berlinsky et al.@ suggested that the Abrikosov lattice at near Hc2 varies continuously from 
a triangular, through an oblique, to a square one with increasing field and s-d mixing 
parameter. In our consideration to the orientation, to smoothly connect the flow trajectories 
of quasiparticles in Fig. ^ with the nearest neighbor vortices, the nearest neighbor direction 
prefers 45°, which is consistent with the result by Won and Maki. 

In the experiment, Keimer et al.0 reported an oblique lattice by a small-angle neutron 
scattering study of the vortex lattice on YBa2Cu307 in a magnetic field region of 0.5T< 
H <5T. The lattice is with an angle of 73° between the two primitive vectors and oriented 
such that the nearest-neighbor direction of vortices makes an angle of 45° from the a-axis. 
The oblique lattice was also observed by scanning tunneling microscopy (STM) by Maggio- 
Aprile et al.ii To understand the origin of the oblique lattice and its orientation, in addition 
to the effect of the intrinsic in-plane anisotropy, that is, the difference of the coherence 
lengths between a-axis and 6-axis directions,^ the study of the vortex lattice in d-wave 
superconductors is certainly important. 

Most probable candidate of experimental probes for detecting the fourfold symmetry of 
the vortex structure is to observe the local density of states (LDOS) by STM. If the four 
ridges placed four-fold symmetrically in LDOS are observed at the core region, it will be 
possible to determine definitely the symmetry of high-Tc superconductors and distinguish 
other symmetries, e.g., the four ridges of the LDOS in ^^.y-wave pairing is rotated by 45° 
from that of the c?a,2_j^2-wave case, a circular ridge of the LDOS is expected around a vortex 
in s-wave pairing. For the energy range 0.05 < E <1, LDOS extends along the x-axis and 
y-axis (or a-axis and 6-axis) as shown in Fig. ^ (a). For E ^ 0, LDOS extends along the 
lines y = ±x as shown in Fig. ^ (d), rotated by 45° from that of higher energy. The peak 
lines are detected as the line regions of the zero-bias conductance enhancement. 
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While the peak hnes for E ^ was previously proposed by Kashiwaya et al.EJ on the 
analogy of bound states in the pseudo-quantum wells and in the surface region of dx2_y2- 
wave superconductors, our explanation of the peak lines is different from theirs. Their theory 
starts from the assumption that quasiparticles rotates around the core along coaxial circles 
as in the case of conventional s-wave superconductors. Along the circular trajectory, all 
quasiparticles feel the node of the pair potential at the point = 7r/4 and its equivalent 
points, where the sign of the factor cos 26' in (i-wave pair potential changes. Thus the 
formation of the mid-gap states are expected at these nodes, and the bound states form two 
orthogonal lines y = ±x. Following their theory, the pair potential is suppressed along the 
lines y = ±x. However, it contradicts with our numerically obtained pair potential which is 
suppressed along the x-axis and y-axis as shown in Fig. |l]. Our understanding is as follows. 
Quasiparticles flow around a vortex along not coaxial circles but four trajectories shown in 
Fig. 1^. Along each of the four trajectories, quasiparticles do not feel the sign change of the 
factor cos2^^. The peak lines y = ±x for E = are due to the quasiparticles approaching 
and leaving the core region with the flow direction 6 = tt/A (and its equivalent directions), 
which is the direction energy gap is suppressed in (ia,2_y2-wave superconductors. We obtain 
the same LDOS also for the case when we use | cos 26*1 instead of cos 26', corresponding to 
the case where the sign of the pair potential does not change. 

The de Haas-van Alphen (dHvA) oscillation in the mixed states is one of the related 
topics with the bound states of quasiparticles around a vortex. The dHvA oscillation was 
reported also on YBa2Cu307.il@ For the dHvA oscillation, enough quasiparticle transfer of 
low energy bound states around a vortex core should occur between vortices. In (ia.2_j^2-wave 
superconductors, the quasiparticles with E = exists along the lines y = ±x. Following 
the theory of Won and Maki,^ the nearest neighbor vortex is located on the lines y = ±x. 
Thus the low energy bound states around each vortex are connected each other by the 
quasiparticle flow with the flow direction (p = tt/A and its equivalent directions. In this 
situation, dHvA oscillation is easier to occur compared with the conventional s-wave case. 
This consideration can be applied to the cases of other symmetry if the energy gap has node. 

16 



Another topic is whether the energy level of quasiparticle bound states is quantized 
around a vortex in d-wave superconductors. In conventional s-wave superconductors, the 
energy level is discretized and the separation of the energy levels is the order of A^/Ep {Ep: 
Fermi energy) .0 Wang and MacDonald^l reported that the separation of energy levels is 
visible in the s-wave case and is not visible in the (i-wave case using a short coherence length 
model. While our quasi-classical calculation can not directly consider the quantization, our 
results of LDOS give the following consideration to this problem. Since quasiparticles do not 
flow along closed circular trajectories in cZ-wave superconductors, the quantizations could be 
different from the s-wave case. For the case of open trajectories as shown in Fig. it is 
expected that the quantization is absent or, if it exists, the separation of the energy levels is 
small, which seems to be consistent with the results by Wang and MacDonald. To consider 
the quantization in detail, we have to consider the coherence effect of quasiparticle wave 
functions at the cross points of four trajectories in Fig. ^. 
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TABLES 

TABLE I. Fitting parameters in Eqs. (3.4) and (3.5) at various temperatures T/Tc = 0.1 ~ 0.7. 
Temperature dependent uniform gap A(r) is also presented. 
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A(r) 


0.1 


0.132 


-3.24 


15.8 


-37.5 


-1.19 


8.74 


-23.7 


1.00 


0.2 
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5.77 
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-4.13 
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-1.38 
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-1.38 
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0.856 


-0.522 


0.648 


-1.26 


-0.0795 


0.285 


-0.717 


0.882 


0.7 


1.03 


-0.468 


0.475 


-1.23 


-0.0492 


0.215 


-0.637 


0.801 
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FIGURES 

FIG. 1. Pair potential at T/Tc = 0.1. It is fourfold symmetric even in the case the mixing of 
the s-wavc component is absent, (a) Contour plot of the amplitude, |A^2_j^2(r)|. From the center, 
0.1, 0.2, • • •, 0.9. Amplitude is suppressed along the 0° direction, (b) Stereographic view of the 
phase, arg Aj.2_j^2(r). It is positive (negative) for < ^ < 7r/4 {tt/A < <p < tt/2). 

FIG. 2. Pair potential as a function of distance from the vortex center, r, at various tempera- 
tures T/Tc = 0.1 ~ 0.7. With decreasing temperature, the size of the vortex core decreases and four- 
fold symmetry becomes clear, (a) Amplitude along the 0° direction, \Ari.2_y2{r,<p = 0)|. (b) Differ- 
ence of the amplitude along the 0° and the 45° directions, \A^2_y2{r, (f) = G)\ — \K^2_y2{r, ^ = 7r/4)|. 
(c) Phase of A^2_y2{Y) along the 22.5° direction, arg A3.2_j^2 (r, = 7r/8). 

FIG. 3. Temperature dependence of fitting parameters ^ (a) and 03, 63 (b). The lines are 
guides for the eye. Core size ^ approaches linearly on lowering temperature. 

FIG. 4. Current distribution around a vortex at T/Tc = 0.1 in units of — 2eiVoUFAo. (a) 
Stereographic view of the amplitude, |J(r)|. There are four small peaks around the core at 45° 
and its equivalent directions, (b) Stereographic view of the radial component, Jr(r). It is negative 
(positive) for < ^ < 7r/4 (7r/4 < ^ < 7r/2). 

FIG. 5. Contour plot of the magnetic field distribution, H{r), at T/Tc = 0.1. The field extends 
along the x-axis and y-axis. 

FIG. 6. Current distribution as a function of r at various temperatures T/Tc = 0.1 ~ 0.7 in 
units of — 2eA^o^^FAo. (a) Amplitude along the 0° direction, |J(r, (/> = 0)|. (b) Difference of the 
amplitude along the 0° and the 45° directions, | J(r, ^ = 0)| — |J(r, = 7r/4)|. (c) Radial component 
along the 22.5° direction, Jr{r, (f) = tt/8). 
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FIG. 7. Magnetic field distribution (arbitrary unit) as a function of r at various temper- 
atures T/Tc = 0.1 ~ 0.7. (a) Deviation from the value of the vortex center Hq along the 0° 
direction, H{r,(j) = 0) — Hq. (b) Difference of the field along the 0° and the 45° directions, 
H{r,(t) = 0) - H{r,(t) = tt/4). 

FIG. 8. Local density of states, N{r,E). Oscillating behavior of peak's height along the peak 
lines is due to the lack of mesh points in these figures. If the number of mesh points is enough 
large, the height smoothly varies along the peak lines, (a) E = 0.2. Self-consistent calculation 
gives the similar fourfold symmetric distribution to that of Schopohl and Maki. (b) The same as 

(a) but wider region is presented. It is noted that the eight peak lines extend outward from the 
vortex, (c) E = 0.01. Peak lines approach the center of vortex, but have the similar structure to 

(b) . (d) E = 0. In addition to the usual large zero bias peak at the vortex center, whose height is 
truncated in this figure, there are four small peak lines on the lines y = ±x. 

FIG. 9. Flow trajectories of quasiparticles with energy E{< 1) around a vortex. Points A, B 
and C correspond to the peaks in Fig. 10. 

FIG. 10. Local density of states as a function of E and r along the direction (j) = vr/S (a), 
^ = (b) and ^ = 7r/4 (c) . Three peak lines A, B and C for £^ < 1 correspond to the trajectory 
in Fig. 9. The peak line at = 1 is the gap edge. Oscillating behavior of peak's height along 
the peak lines is due to the lack of mesh points in these figures. If the number of mesh points is 
enough large, the height smoothly varies along the peak lines. 
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